Nematic-Isotropic Transition with Quenched Disorder 
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Nematic elastomers do not show the discontinuous, first-order, phase transition that the Landau- 
De Gennes mean field theory predicts for a quadrupolar ordering in 3D. We attribute this behavior 
to the presence of network crosslinks, which act as sources of quenched orientational disorder. We 
show that the addition of weak random anisotropy results in a singular renormalization of the 
Landau-De Gennes expression, adding an energy term proportional to the inverse quartic power of 
order parameter Q. This reduces the first-order discontinuity in Q. For sufficiently high disorder 
strength the jump disappears altogether and the phase transition becomes continuous, in some 
ways resembling the supercritical transitions in external field. 
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I. INTRODUCTION 

Although there is a large volume of literature devoted 
to the effects of quenched disorder, there has been rel- 
atively little study on how it influences the behavior of 
systems whose pure versions undergo a first-order phase 
transition. This question was first addressed by Imry and 
Wortis who showed that inhomogeneities may cause 
local variations of the transition temperature inside the 
sample. Provided that the cost in interface energy is not 
great, bubbles of the "wrong" phase are formed, eventu- 
ally leading to a substantial rounding of the transition. 
A theorem due to Aizenman and Wehr |2| shows that in 
less than two dimensions there can be no phase coexis- 
tence at the transition, and therefore no latent heat, in 
a system with quenched random impurities. Therefore 
these systems are expected to always exhibit a continu- 
ous transition. 

The influence of quenched impurities coupling to the 
local energy density has been extensively studied by 
Cardy 0,0- A mapping to the random field Ising model, 
whose renormalization group flows are well known, ad- 
dresses the question of what happens in higher than two 
dimensions where the Aizenman- Wehr theorem is not ap- 
plicable. It was found that, depending on the specific 
values of parameters such as the laten heat and the sur- 
face tension, the phase transition can either be first or 
second order. More conclusive analytic results were ob- 
tained in an Ising model with discrete order parameter 
0. The pure system exhibits a fluctuation-driven first 
order transition, where the mean field theory predicts a 
continuous transition but fluctuation effects make it dis- 
continuous. Quenched randomness eventually drives the 
transition to become continuous in two dimensions, in 
accordance with the Aizenman- Wehr theorem, but this 
may or may not happen in higher dimensions. 

The majority of studies in this field are carried out for 
spin-glass or analogous systems, e.g. with a frustrated 
dipolar ordering However, in such systems the ex- 



perimental work is difficult and results are often indirect. 
In contrast, the quadrupolar orientational ordering of ne- 
matic liquid crystals offers the easy experimental access 
to thermodynamic and structural features of phase tran- 
sitions. The classical work on frustrated nematic liquids 
by Bellini et al. 0, 111 has generated a large interest in 
studies of liquid crystals in random environments, such 
as porous silica gel. One must appreciate, however, that 
the characteristic length scale of such disorder is much 
greater than the coherence length of the nematic order 
parameter and thus the theoretical concept of a continu- 
ous coarse-grained random field 9] is difficult to sustain. 

Quenched disorder is intrinsically present in nematic 
elastomers as a direct result of their synthesis 0, ^] . 
Sources of disorder are introduced by crosslinking a 
liquid-crystalline polymer melt. In the simplest situa- 
tion, the crosslinking takes place in the isotropic phase, 
in which case the local anisotropy axis of each crosslink- 
ing moiety is randomly oriented. Once the polymer net- 
work is formed, the configuration of the crosslinks re- 
mains quenched, that is, it does not change with time and 
temperature. Unless special precautions are taken dur- 
ing network fabrication, the low temperature (ordered) 
phase of nematic elastomers is always an equilibrium 
polydomain director texture [T^ , . This is in marked 
contrast with a kinetic "polydomain" texture often re- 
ferred to as Schlieren texture 01 > which is the conse- 
quence of nucleation and growth mismatch in a system 
undergoing the first-order transition |15|. The equilib- 
rium polydomain structure of nematic elastomers is re- 
versible with changing temperature and is characterized 
by the uniform non-zero order parameter, but the highly 
non-uniform orientation of the principal axis of nematic 
director n. Correlations between directors decay rapidly 
and eventually vanish at distances much larger than £, 
the correlation length or domain size. This is in agree- 
ment with the general result that quenched impurities 
destroy long-range order, first shown by Larkin |l6j and 
then generalized by Imry and Ma Q . There is a full anal- 



2 



ogy with a corresponding dipolar system named "random 
anisotropy magnets" 0,0,0- In fact, all other (e.g. 
smectic |2I|) liquid crystal elastomers follow the same 
pattern of forming the equilibrium textures with a char- 
acteristic length scale often referred to as the domain 
size. 

This length scale in typical nematic elastomers is of the 
order of microns 0, E| , therefore light passing through 
the sample is multiply scattered on birefringent domains 
with randomly oriented optical axis [2l| (see |Hj for a 
brief review of experimental facts in this area). As a 
result such a sample is completely transparent at high 
temperatures, but becomes opaque below its nematic-to- 
isotropic transition temperature Tjvj. It should be noted 
that the polydomain texture is the thermodynamically 
stable low-temperature phase. Applying an adequately 
strong aligning stress |l2( or increasing the temperature 
above Tni destroys the polydomain texture, by align- 
ing the local axis of each domain, or removing the opti- 
cal contrast between them. However, once the stress is 
removed (or the temperature lowered), the elastomer re- 
turns to its previous state and the average domain size is 
found to be reversible during this stress (or temperature) 
cycling. 

Experimental investigations in nematic elastomers 
support the theoretical ideas of a continuous nematic- 
isotropic transition, rather than the discontinuous first- 
order transition expected by the Landau-De Gcnncs 
mean field theory based on the symmetry of quadrupolar 
ordering in three dimensions. The fact that polydomain 
nematic textures are optically opaque creates a practical 
problem when attempting to experimentally determine 
the local order parameter Q of the mesogenic units. It is 
impossible to use birefringence, dichroism or Xray mea- 
surements, the methods that have made this task simple 
in aligned liquid-crystal systems. Nuclear magnetic res- 
onance (NMR) provides perhaps the only opportunity of 
order detection, by tracking the bias in orientational mo- 
tion of selected chemical bonds and providing a unique 
Q-signature even in the orientationally averaged case. It 
was used to measure Q(T) of both a nematic polymer 
melt and its corresponding crosslinks network p2l l23l| . 
The addition of crosslinks was shown to make the nematic 
transition smooth, as well as to slightly reduce order. An 
indirect alternative to NMR comes from applying an ex- 
ternal field to align the domains. Birefringence can then 
be measured, as long as the system has passed the criti- 
cal point of the polydomain- monodomain transition . 
The order parameter can thus be found for a series of 
decreasing applied fields. Although the zero-stress limit 
is not accessible, it can be extrapolated from the compar- 
ison of the other curves and, again, a continuous phase 
transition is seen j24[. 

The continuous transition is also found in carefully syn- 
thesized monodomain nematic elastomers, where a sec- 
ond stage of the crosslinking takes place when the sample 
is stressed 25] . Plots of the macroscopic order parameter 
Q(T), obtained through birefringence [2(|, X-ray scat- 



tering and NMR measurements [2^, show the same 
behavior. There are two possible explanations for such a 
deviation from the basic symmetry-based expectation of 
a first-order transition. Inhomogeneities may cause local 
variations of the transition temperature inside the sample 
and lead to a substantial rounding of the transition, as 
discussed in . The other explanation considers stresses 
imprinted in the system during the second crosslinking 
stage, which add a —fQ term to the Landau-De Gennes 
expansion. Obviously, for adequately large /, the system 
would become supercritical and show a continuous tran- 
sition. The analysis of the NMR spectra supports the 
supercritical scenario, although some degree of inhomo- 
geneity was also found using the second moment of the 
spectra [23] ■ Another study of strain as a function of tem- 
perature over a range of applied tensile stress argues the 
opposite point [2^ . The blurring of the isotropic-nematic 
transition is attributed to the presence of heterogeneities, 
although the authors do not consider boundary effects 
that are bound to be present as discussed in p|. 

Uchida [2^] has studied disordered polydomain ne- 
matic elastomers with emphasis on the role of nonlo- 
cal elastic interactions. He has shown that networks 
crosslinked in the isotropic phase lose their long-range 
orientational order due to the locally quenched ran- 
dom stresses, which were incorporated into the affine- 
deformation model of nematic rubber elasticity. Sim- 
ulation work was carried out to investigate the role of 
random bonds and random fields that might be present 
in elastomers |30|, in both cases finding that the first- 
order isotropic-nematic transition to broaden into a 
smooth crossover. For random-field disorder, the smooth 
crossover into an ordered state is also attributed to the 
long-range elastic interaction present in elastomers. A 
recent coarse-grained model for liquid-crystalline elas- 
tomers has also found that both homogeneous and in- 
homogeneous samples undergo a continuous isotropic- 
nematic transition |3l| . 

In this paper we apply the traditional spin-glass tech- 
niques to investigate the characteristics of nematic phase 
transition in a system with quenched random anisotropy. 
The structure of the paper is as follows. In Section II 
we summarize a physical model of quenched disorder in 
nematic systems following |l9t , and introduce the replica 
Hamiltonian. Section III applies the auxiliary fields to 
incorporate several constraints into this problem and ob- 
tains the effective mean-field free energy of disorder in 
the system. In Section IV we investigate the stability 
of replica symmetry and discover the limits where our 
solutions are valid. Finally, in Section V we obtain the 
final free energy renormalization in terms of the order 
parameter Q and disorder strength, and investigate the 
characteristics of the nematic phase transition in various 
situations. We conclude by discussing our results and 
comparing them with experiments. 
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FIG. 1: Schematic representation of how crosslinks provide 
easy anisotropy axes {k}. The nematic director is forced to be 
aligned, in the vicinity of the crosslink, with the axes, which 
are represent by the arrows. Both the orientation of {k} as 
well as the positions of the crosslinks {Ri} are random. Since 
the crosslinks are confined by the network topology, they add 
quenched disorder to the nematic system. 



II. MODEL 
A. Sources of Quenched Disorder 

In the case of nematic elastomers crosslinked in the 
isotropic phase, the sources of quenched disorder are pro- 
vided by the network crosslinks. Almost independently 
of their specific chemical structure, the crosslinks contain 
anisotropic groups that locally provide easy anisotropy 
axes k: it is favorable for the local director to align along 
k in the vicinity of the crosslink because the anisotropic 
molecules in both the crosslinkers and in the nematic 
system interact, both sterically and via the long-range 
van de Waals attraction. The local anisotropy axes of 
the crosslinks, as well as their distribution inside the 
sample, are quenched variables since the crosslinks can 
neither rotate nor move once the chemical synthesis of 
the elastomer is complete. Although this has never been 
tested experimentally, there are two independent molecu- 
lar models that estimate the energy of orientational con- 
finement that a crosslink experiences from the surround- 
ing network strands |32t l33j 

We follow earlier work [l(| in modelling the local cou- 
pling of the nematic order and the random field applied 
by the crosslinks. For a crosslink positioned at Ri, with 
an anisotropy axis ki, an energy — ^yki ■ Q ■ fcj is raised 

due to the interaction with the local nematic order pa- 
rameter Qij — Q (jiiUj — 6ij/3), where 7 is the coupling 
strength. Employing a coarse-grained expression for the 
continuum density of crosslinks p (r) = 5(r — Ri) 

and substituting the full tensor expression for Q, we get 
the random field energy: 



J Sr 1 Qp{r) (k-n) 2 , 



(1) 



where the irrelevant constant 7 J p(r) has been dropped. 
Although the random field energy is proposed specifically 
for a nematic elastomer, it is a general expression which is 



valid for all random-anisotropy systems with underlying 
quadrupolar symmetry. 

To obtain the full Hamiltonian describing the nematic 
ordering, the gradient elasticity penalizing the fluctua- 
tions of the director field must be also taken into account. 
In continuum elasticity the Hamiltonian takes the form: 



H\p,k] = Jd 3 



y (Vn) 2 - 1 Qp{k-nf 



(2) 



where K is the Frank elastic constant in the one-constant 
approximation. A simple dimensional argument gives 
K ~ ksT/a, where a is the nematic coherence length, 
below which the meanings of the director n and or- 
der parameter Q are ill-defined |l4j . Both microscopic 
and phenomenological theories of nematic-isotropic tran- 
sition give the elastic constant K to scale as Q 2 for 
Q <C 1. Combining these two estimates, we take that 
for small Q the elastic constant is approximately given 
by K ~ kBTQ 2 /a. It is important to clarify that we 
are examining a homogeneous sample and as a result the 
magnitude of Q and the Tmi are uniform across the sam- 
ple. There is a rich literature on the role inhomogeneities 
play in general first-order systems [l|| and more specifi- 
cally in nematic elastomers [5^, l29i I30L l34l |. The im- 
portant difference in our assumptions is that, although 
the director correlations in the polydomain nematic are 
short-range (equilibrium spin- glass texture), the under- 
lying nematic order parameter Q is homogeneous across 
the system. This assumption is based on the fact that the 
spin-glass like nematic textures are in fact homogeneous, 
in the sense that every element of the sample is equivalent 
to others: no 'domain walls' (unlike, for instance, during 
the stress-induced poly-monodomain transition |10| when 
the domain walls localize). 



B. The Replica Method 

There are three established methods of dealing with 
quenched random fields in the replica method framework: 
one based on the functional renormalization group anal- 
ysis, another using the Gaussian variational method and 
the third using auxiliary fields. This paper employs the 
latter. 

We are interested in results that do not depend on 
the specific distributions of p(r) and k because we can- 
not control these distributions experimentally. In other 
words we are looking for the free energy averaged over 
the random distributions of the quenched variables p(r) 
and k. Crosslinking the sample above Tjvj makes the 
easy anisotropy axes point at random directions, with 
an isotropic probability of orientation P(k) = 3—. Fur- 
thermore, the crosslinks are dispersed randomly in the 
sample with density p. The probability that a particular 
distribution p(r) occurs is Gaussian: 



P[p(r)] 



exp 



2po 



(3) 
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where po is the mean density of crosslinks. 

It is not possible to directly average the logarithm of 
the partition function Z and obtain the exact free en- 
ergy. So an alternative definition of a logarithm (the 
limit: logZ = d m Z m \ m ^o) is used, allowing to per- 
form the simpler average of the product Z m . This way 
of dealing with quenched disorder is called the "replica 
trick" , first introduced in the context of spin glasses by 
Edwards and Anderson (35|. The expression for the free 
energy arising from disorder then reads: 



F d = -k B T (log Z) pM = -k B T^-(Z 



dm 



(4) 



= -k B T 



_d_ 

dm 



m „ 

exp [-ftH Iep ] 

»— a=1 J 



where we now have m identical "replicas" of the system, 
labelled by the index a. The aim of this work is to ob- 
tain Fd as a function of the order parameter Q and add 
it to the Landau-De Gennes free energy to see how it 
influences the phase transition. 

A rough sketch of the averaging over disorder is given 
below. The density average over P(p) yields a random 
field term 



exp 



(/3jQ) 2 p (k a ■ n a ) 2 (k b ■ n b ) 



where (here and throughout this paper) (3 = l/k B T. The 
reader's attention is drawn to the appearance of a second 
replica index b due to the square of the J2 a - Since the 
distribution of the orientations of the easy axes P{k) is 
assumed fully isotropic, symmetry arguments show that 
the average (kikjkik m )p(k) is proportional to (6ijSi m + 
duSjm +dim.3ji)- The constant of proportionality is equal 
to l/(2d + d 2 ), where d is the dimensionality of the unit 
vector k. Therefore 



( e )P{k) = 

= 1 + (Wk l k m ) rin 3 n l n m + 

~ „(k z & k L k m ) n z n j n l n m 



(5) 



Since (n a - rib) < 1 higher terms of the Taylor expan- 
sion are smaller than the lowest order term ({n a ■ n b ) 2 ) 
and are subsequently dropped in the last line of Eq. JSJ. 
This approximate treatment retains a random field term 
which is overall fourth order in n, as it was before the 
fc-averaging, and it is most frequently met, and used, in 
molecular theories involving rotational diffusion. 

The second line of Eq. (0J is obtained after averaging 
over quenched disorder, sketched above, and provides the 
definition of the "Replica Hamiltonian" , which no longer 
depends on the quenched distributions of {k} and p(r), 



but instead couples different replicas of the system: 

m r ( k 
H rcp [n (r)] =J2 d3r \ Y ( V " Q ) 2 5ab + 



2 (n a ■ n b ) + (n„ • n a ) (n b ■ n b ) 



(6) 



where subscripts a and b are the replica indexes and m 
is the number of replicas that will be set to zero at the 
end of the calculation. Parameter T, arising from com- 
pleting the Gaussian square between the Eq. © and the 
random-field term in the Eq. (J2J, reflects the strength 
of the disorder and has a quadratic dependance on the 
order parameter: 



r 



7 2 Po 
15 k B T 



(7) 



It is noted that all replicas are assumed to have equal 
disorder strength and equal magnitude of the local order 
parameter, i.e. j a = j b and Q a — Q b for all a and b. 
Bearing in mind that the director is a unit vector we see 
that the term with (n a ■ n a ){n b ■ n b ) in the random field 
part of the Replica Hamiltonian just adds an irrelevant 
constant to the expression. We drop this term and keep 
the relevant contribution — T (n a ■ n b ) 2 that describes the 
coupling between different replicas. 



III. DISORDER FREE ENERGY 

A. Auxiliary Fields 

Care must be taken to ensure that the director n(r) 
remains a unit vector. Although this was assumed to be 
the case, it has not been implemented explicitly in the 
Eq. 10. One way to achieve this multiplies the partition 
function in Eq. with the delta-function constraint 



6(n 2 



1 



1 

27 



(8) 



where <p is an auxiliary field that allows the delta-function 
to be written in its exponential form. We proceed with a 
mean-field treatment of the auxiliary fields where a con- 
stant value for </> is assumed, independent of spacial coor- 
dinates. This approximation implies that the constraint 
n 2 = 1 is equally enforced across the whole sample. It is 
reasonable to expect this, given that the sample is spa- 
tially homogeneous. The same reasoning explains why <f) 
has no dependance on the replica indexes: since its value 
is independent of the position of the crosslinks in the 
sample, it cannot have different values for different repli- 
cas. The corresponding quadratic term J2 a * ( n a ~ -0 
has to be added to the Replica Hamiltonian in Eq. 10. 

To obtain the disorder energy one must evaluate the 
statistical sum over all possible trajectories n a . The stan- 
dard way to evaluate Hamiltonians with quartic interac- 
tions is to introduce an auxiliary field, here a tensor X ab , 
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which reduces the Hamiltonian to bilinear order in n a . 
To employ the method all the quantities in /3H Tep [n(r)] 
must be dimensionless. Since the integral over r has the 
dimensions of volume, we move to a discrete summation 
over all points in space: 



dx 



dy 



dz 



E ■ 

points r 



The limits of the r-integration are L, the size of the sys- 
tem, and the short-distance cutoff a - the nematic coher- 
ence length, below which the continuum representation 
is no longer applicable. The (a, b) replica coupling term 
then becomes 



exp 



f3Ta 3 ( n a ■ n b ) 



(9) 



J dXab exp < ^ 



+ Aafc(«a ' Ub) 

4r 



which involves the dimensionless constant f = (3Ta 3 . It 
is important to clarify the meaning of the A 2 ^ term in 
Eq. 10: it is the square of the value X a b rather than 
an element of the product of two matrices. Further- 
more, from now on we shall use a mean-field approxima- 
tion, where it is assumed that X a b has no r-dependancc. 
This is an important limitation, but we believe it is rea- 
sonable as we are looking for homogeneous ordering in 
the system. Summation over r of the \ 2 ah term yields 
TV A 2 b /(4r), where N = V/a 3 is the number of "discrete 
spacial points" and V = L 3 is the system's volume. 

Moving to the corresponding discrete Fourier space, 
the effective Replica Hamiltonian includes both auxiliary 
fields: 



PH eB [n (q)] = 



icj)5 al 



A 2 

+ N-^ + 

4r 



E 



9 L 



Kq 2 



+ i4> Sab - X a t 



(n„-n 6 ) }, (10) 



with the dimensionless elasticity constant K — (3Ka 3 . 
The discrete sum over q is related to the integral via 
= L 3 J p^yg ■ As mentioned above, the conversion 
from the integral to the discrete sum is essential so that 
all the quantities in H c g (such as T, <fi, X ao and Kq 2 ) 
remain dimensionless and the logarithm of their sum can 
be evaluated correctly. 

To be able to deduce the disorder energy as a function 
of the order parameter, we make explicit the dependence 
of r and K on the magnitude of the order parameter Q: 



f = gQ 2 , with 
K = Pna 3 Q 2 , 



(7/3)W 
15 



(11) 



where k ~ fe^T/a is the Frank elastic constant deep in- 
side the nematic phase. It is worth noting that T is always 



significantly smaller than one. The distance between 
crosslinks, d c , which can be deduced from the crosslink 
density po — d~ 3 , is found around 7 — lOnm in nematic 
elastomers 10] . A typical coherence length for nematics 
is a ~ 5nm, hence (a/d c ) 3 ~ 0.4 or less. The coupling of 
disorder to the nematic director is deduced from the size 
of the domains in Ref ^(| and is found to be 7 ~ OAUbT. 
Therefore this crude estimate gives g ~ 4 x 10~ 3 for the 
nematic elastomers studied in the literature. 



B. Replica Symmetry Case 

To be able to advance in the calculation, a matrix form 
of the auxiliary field A a b has to be postulated. A reason- 
able starting point is to assume that it has the simplest 
possible form, where its elements have a constant value 
independent of a and b: 



Xab — X (tab ~ 5 a b) ■ 



(12) 



where 1„& is the matrix with all its elements equal to 
one and S a b is the identity matrix. This scheme is fre- 
quently encountered in the literature [3jj and is called 
the "Replica Symmetry" limit. It is important to clar- 
ify that we are free to choose any form we like for X a t, 
and that we will later come back to this choice and check 
whether it is appropriate or not, in Section llVI There 
is an important reason why the diagonal elements of the 
auxiliary field A a f, are chosen to be zero in Eq. |Q . Had 
X a b = Xt a b been our choice, then the summation over a 
and b of X 2 b would have given m 2 X. The quadratic depen- 
dence on m would have meant that, after differentiating 
with respect to m and setting m = 0, this term would 
be equal to zero. Clearly this is not acceptable since the 
introduction of the auxiliary field in Eq. JHJ) requires a 
non-zero X 2 ab quadratic term. 

Substituting the replica-symmetric ansatz into the ef- 
fective Hamiltonian of Eq. we find: 



0H eS = J2 { 



A 

-icp6ab + N^(l ab -S ab ) + 

4r 

4 E^K 



(13) 



where the propagator of the a — b replica coupling is given 
by: 



g: 



Kq 2 + 2i(f) + 2A 5ab - 2A1 



(14) 



A consequence of replica symmetry is that only in- 
volves matrices 5 a b and l a &, and as a result its logarithm 
is easily obtained: 



lab (15) 



logG-, 1 = log [Kq' + 2icf> + 2X)5 ab + 



1 , / 2Am 
— log 1 - -~— 

to \ Kq 2 



2iq 



2X, 
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The path integral over configurations n a with the 
statistical weight determined by the effective Hamil- 
tonian H13JI is Gaussian and gives (DetG -1 ) -1 / 2 = 
cxp(— A tr log G _1 ) for each of the three vector compo- 
nents of n. As a result the disorder free energy is given 
by: 

m = --— Jdxjdd, exp | - -^(to 2 - m) + 



+id>m + 



(16) 



from the three (identical) path integrals for the compo- 
nents of n. 



C. Disorder Free Energy 

The aim of this section is to determine the particular 
values of the auxiliary fields (A* b and dj*) that make the 
disorder energy of Eq. i|16|) a minimum. To treat the 
problem properly, one would have to evaluate the inte- 
grals over X a b and </>, which is analytically challenging. 
The standard way to bypass this difficulty, is to employ 
the saddle-point approximation based on the simple ob- 
servation that the exponentially most significant contri- 
bution in Eq. I|16|) will occur when the exponent is a 
maximum. Therefore 



pF d (K, T) « -i^-exp<{ min 

2 dm [ \,ic/> 



i5>logG a t 



+im(j) - 



iVA 



4r 



— (to — to) 



(17) 



where min,\,i0(..) represents the minimum of a function 
with respect to variations in A and After substituting 
the trace of the logarithm from the Eq. Q15[l. we differ- 
entiate with respect to to and then set to equal to zero. 
The disorder energy then takes the form: 



j3Fd ss min ^ — ic, 

A, i(f> 



-E 

2 ^ 

q 



log (Rq 2 



NX 2 
4f 

2id> + 2X 



(18) 



2A 



Kq 2 + 2i<j) + 2X 



and we are left with the task of finding the stationary 
point (idj* , A* ) . From now on we move to the contin- 
uum limit of space, where the discrete sum is replaced by 
J2 q ~~ * V J Hnjz d<l an d the coherence length a is taken 
to zero limit. 



1. Optimal i<f> in the absence of disorder 

The auxiliary field </> ensures that the nematic direc- 
tor is a unit vector. This constraint should of course be 



enforced whether the disorder is present or not. In fact, 
there is no physical reason why the inclusion of disorder 
should significantly alter this constraint. Hence, as a first 
approximation, wc look for the optimum idj* when there 
is no disorder in the system. Setting both T and A equal 
zero in the Eq. IjlSJ) and differentiating with respect to 
id> we obtain: 



didj 



= 



(19) 



1 = 



2vr 2 



1^4 _! 




In the continuous limit of (7, nax — » oo and the arctangent 
is equal to 7r/2. This equation can be re-written as 

(20) 



-1 



2tt 2 K 



Clearly the —1 term is negligible compared to 
(QmaxV/K), and can be therefore neglected. Another 
factor that supports this omission is that we want to ex- 
amine what happens close to the phase transition, where 
Q — * 0. In this limit, the elastic constant is known to 
vanish (oc Q 2 ) and the term with K^ 1 dominates in the 
bracket. 

Assuming that the inclusion of disorder has only a mi- 
nor effect on d>, the value we will use from now on is: 



i(f>* = 



2Kq 



2 

max 



(21) 



which tends to infinity in the continuum limit of space. 
It is interesting to note that very close to the transition 
the constraint relaxes since dj — > 0. This is not surprising 
since the meaning of the director itself becomes ill-defined 
as we approach the transition point. 

The full calculation to obtain <^>*(A) in a system with 
disorder is possible. However, using a disorder-dependent 
4>*(X) in Eq. (|18|) makes it analytically impossible to solve 
dFd/dX — and determine the optimum A* . Even a per- 
turbative approach : dj* (A) = <^ =0 + " S77la ^ correction" 
does not help matters. To overcome this difficulty the 
(</>*, A*) saddle point should be found numerically, search- 
ing for a global energy minimum in the Q — 7 space. The 
important drawback of the latter is that we will then be 
unable to determine the analytical form of X*(Q) and 
therefore the final Fd(Q) correction to the Landau free 
energy of the phase transition. 



2. Optimal A for weak disorder 

We proceed to determine the value of the auxiliary 
field A* that minimizes the energy in the replica symmet- 
ric approximation. Differentiating the right-hand side of 
Eq. I|18|) and demanding it to be zero we find: 

2 dq. 



X*N 



2T 



3V 



(Kq 2 



+ 2i<j>* + 2X" 
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Integration over momentum space gives the stationary 
condition on the auxiliary field: 

A*iV 3V 



2T 




(2K)V 2 v^+I 



In the continuous limit of g max = 2ir/a — > oo the first 
term vanishes and the arctangent is equal to tt/2. Equa- 
tion 1(221) has only one real solution. Unfortunately, its 
full expression is too long and cumbersome to appear 
here explicitly; instead we demonstrate its behavior in 
two limits. Expanded in powers of T < 1 (weak disor- 
der) it takes the form 

A*(r) = — _ — r + o(r 2 ). (23) 



It is a reassuring property that A* vanishes as T — ► 0. 

Critical to our work is the behavior of A* as the order 
parameter Q tends to zero. Both T (from its definition) 
and K (for small Q) are quadratic functions of Q. There- 
fore the leading term in the series expansion of small Q 
is 



A*(Q) 



(3V) 2 / 3 f 2 / 3 



27r 2/ 3i V2/3 K 



0(Q 4/3 ). 



(24) 



The scaling A* cx Q~ 2 / 3 is thus obtained, showing that 
A* diverges as the transition is approached, that is, even 
a weak disorder becomes relevant near the Q — ► point. 



3. Final disorder free energy 

To find the final disorder energy the values of fields 
cf>* and A* are put back in the Eq. (|18J) . Performing the 
q-integrations in the continuum limit we obtain: 



F d 



NX* 
-id> — 

4r 

3V2tt (A 2 



V 

+ 12^2 

idj*\* - 



6g ma x(A* -2i</>*) 

k 



2d>* 



Ky/id)* + X* 



(25) 



The energy we are interested in arises from disorder and 
we can safely ignore terms that are still present when 
A* = r = 0. At diminishing order parameter, the leading 
term of Eq. (|25|1 takes the form: 



F, = 



V(A*) 3 / 2 3V 2 f 
2V2~irK 3 / 2 ~ 8tt 2 K 3 N 



oc Q~ 



(26) 



which clearly diverges as Q — > 0. This divergence of the 
disorder free energy implies that the isotropic phase can 
never be reached and, as we shall see in greater detail 
in Section This in turn leads to the rounding of the 
nematic-isotropic phase transition. Before all this is dis- 
cussed, we examine if the replica-symmetric form of X a b 
in Eq. l|12f) was an appropriate choice, which is far from 
obvious. 



IV. STABILITY OF REPLICA SYMMETRY 
A. The Hessian 

A necessary condition for the replica-symmetric solu- 
tion to be applicable is that the disorder energy is stable 
for infinitesimal variations of that solution. The aux- 
iliary field A* gives an energy extremum; whether this 
extremum is a maximum or a minimum is determined 
by a stability analysis following the work of de Almeida 
and Thouless in spin glasses [37J. We start by allowing 
the matrix of the auxiliary field X a b to deviate from its 
replica-symmetric form: 



A a fc = A* (tab — Sab) + tab 



(27) 



with tab the arbitrary infinitesimal deviation from replica 
symmetry. The disorder free energy will be expanded to 
second order in e: 



Fd(\ a b) — FrS + r ^ H a d,bc e ab£cd 



(28) 



abed 



where Frs represents the energy of the replica-symmetric 
case (e a b — 0), given by the Eq. (|25(l . The second-order 
term in e a b involves the fourth rank tensor of coefficients 
Hadbci which is called the Hessian. It plays an important 
role in this analysis, because the replica-symmetric solu- 
tion is only stable as long as H a dbc is positive definite, 
or equivalently only as long as its eigenvalues are non- 
negative. Of course there is no term linear in e a b because 
its coefficient is given by the right hand side of Eq. 122(1 
and is therefore equal to zero. 

We proceed to find the Hessian of this model. Expand- 
ing the first term of the disorder energy in powers of e a b 
gives: 



N 2 _ N 

— ~^ab — 

4r 4r 



r ( A ° 



rs 2 
b 



+ 2A^e a fe + Cabf-ab 



(29) 



where A^ = A* (l a b — 5 a b)- Hence the contribution to 
the Hessian from this term is (N/2T) S a d 5b c - 

The other contribution comes from the trace of the 
logarithm. The propagator of Eq. ((14(1 can be written as 

G^ 1 = i^Kq 2 + 2id) \ S a b~ 2X a b, therefore allowing A a & to 
vary gives the following form to the propagator 



g: 



Gr 
n 



2e. 



ah ) 



(30) 



where the replica-symmetric G r a s b 
Therefore, 



is given in Eq. 114(1 . 



tr log G-J = tr log [g^ 1 (6 cb - 2G r cd e db )] (31) 

where summation over the dummy indexes c, d is implicit 
and 



Grs 
ab — ~ 



1 



2A* 



Kq 2 + 2i(j)* + 2X* 



(Kq 2 + 2i4>* + 2X*) 2 



"■ab 
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is the inverse of Eq. I|14|) . In the case of two commuting 
matrices: \og{A • B) — log(A) 4- log(-B). It turns out 
that the eigenvectors of the Hessian (e a b) indeed commute 
with matrix G a b- This is largely because the latter is a 
combination of two very simple matrices l j and 6 a b- 
Breaking up the product under the logarithm: 



logG- 1 = logG r 



log(<5 - 2 G rs • e) 



and expanding the second term in a Taylor series, we 
obtain: 



tr log G-J = tr log G a l~ 1 - 2 £ G r b s a e ab 



abed 



(32) 



There are three identical traces to be considered, one 
for each component of the nematic director. Taking into 
account the 1/2 factor in front of the Hessian, this con- 
tribution is —6 J2 q Gda and the overall Hessian takes 
the form: 



N 

^~ SadSbc - 6 G r da G 



(33) 



This Hessian has three distinct values of its many ma- 
trix elements, depending on how its many indexes are 
common to the pairs {a,d} and {b, c}. These values are 

H aa ,bb = P, Haa.bc = H a d,bb = Q and H a( i he = R- 



B. Hessian Eigenvalues 

The eigenvalues of such a 4-rank tensor have been com- 
puted by de Almeida and Thouless in their classical work 
on replica symmetry breaking in spin glasses |37j |. For 
m — > only two distinct eigenvalues exist: 

Ai = P - 4Q + 3R and A 3 = P - 2Q + R (34) 

and it is straightforward to obtain them by inserting the 
appropriate forms of G rs into the expression of the Hes- 
sian of Eq. I|33|) . The first eigenvalue 



Ai = 



2T J (2ir) 3 



1 



{Kq 2 + 2i(f> + 2A) 2 
4A 



N 



(Kq 2 + 2i<t> + 2A) 3 
3(2i(f> + X)V 



dq 



2T 8\/2~7rK 3 / 2 (i(f) + A) 3 / 2 



(35) 



is degenerate and corresponds to two eigenvectors e\ and 
€■2- The first, Ci, is symmetric under interchange of in- 
dexes (e Q b = a for all a, b) and determines whether the 
replica-symmetric fixed point of Eq. H22|) is stable or not. 
In other words, if Aj < the replica-symmetric solution 
corresponds to an energy maximum and has no physical 
relevance. 



The remaining two eigenvectors check the general sta- 
bility of the replica-symmetric scheme. Contrary to e\, 
the second eigenvector e2, corresponding to the degen- 
erate Ai, is symmetric under interchange of all but one 
index (e a b = j3 for a or b = c and e ab = 7 otherwise). 
This eigenvector is not symmetric and as a result a nega- 
tive Ai also means that replica symmetry must be broken 
to determine the correct A a h. The second eigenvalue is 
non-degenerate and is given by: 



A, 



N 
2T 

N 



6V 



A-Kq 2 



1 



(2tt) 3 (Kq 2 + 2icf> + 2\) 2 
W 



dq 



2f 4\72ttA' 3 / 2 (^ + A)V2 



(36) 



and its corresponding eigenvector, called £3, also does 
not have a symmetric form (e a h = S for a = c and b = d, 
( ab = C f° r a = c or « = rf and b ^ c, d, e a b = V otherwise). 
Similarly to the previous case, if A3 < then the replica 
symmetric solution breaks down and different forms of 
X a b must be sought. 

Accordingly, a sufficient condition for the replica- 
symmetric solution to be stable is that all the eigenvalues 
of the Hessian remain positive. Both Ai and A3 show 
a similar behavior as functions of disorder strength 7, 
which was defined in the text above Eq. QJ, and the or- 
der parameter Q. Substituting the optimal (jf and A*, we 
find that when Q does not tend to zero, both eigenvalues 
remain positive, especially for weak disorder, r< 1. As 
the continuous phase transition is approached and the or- 
der parameter diminishes, the three parameters appear- 
ing in Eq. Ij3f)|l become simple functions of Q: r a Q 2 , 
K oc Q 2 and A* esc Q -2 / 3 . Hence, the eigenvalues scale 
as Q~ 2 — Q~ 8 / 9 and, therefore, become negative at some 
point Q 1. This means that for Q — > the replica- 
symmetric solution of the auxiliary field Eq. I|24|) and 
the resulting expression for the disorder energy, Eq. 12 (it . 
correspond to an energy maximum and should not be 
used. 

To find the exact point at which replica symmetry 
becomes unstable we note that it is the second (non- 
degenerate) eigenvalue, A 3 that becomes negative first as 
Q — * 0. Substituting the parameters T, K and A* we find 
this crossover value for stability: 



/stab 



4\/27r(/3Ka) 3 / 2 



(37) 



where the relation V/N = a 3 has been used, with a the 
short-length cutoff. Bearing in mind that g oc a 3 and that 
rea 1/a we find Q s tab oc a 3 . To be consistent with previ- 
ous calculations the continuum limit of space is employed 
and a — * 0. Therefore the threshold value below which 
replica symmetry breaks tends to zero and the Eq. (22} 
is always valid in the continuum limit. Another way to 
see this is by examining the explicit dependance of A3 on 
the cut-off length-scale a. The first term (N/T) scales as 



whereas the second term scales as 



Therefore 



when a is taken to zero the first term dominates and the 
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eigenvalue remains positive. To complete the stability 
analysis, we note that the other eigenvalue also becomes 
negative for smaller values of Q when Q < Q s tab/2, but 
this is irrelevant since the instability of A3 occurs first. 

V. PHASE TRANSITION ANALYSIS 

The aim of this paper is to discuss how quenched orien- 
tational disorder affects the phase transition of nematic 
systems. The total free energy density is a combination of 
the Landau-Dc Gennes expansion of the order parameter 
plus the disorder part (Fd/V) obtained in Eq. To 
get the total energy as an expansion of the order param- 
eter only the leading order contribution of Fd in Eq. i|26|) 
is considered: 

F =^(T-T*)Q 2 -jQ 3 + ^Q A + jQ-\ (38) 
where 

D = 2^r kBT m 

and A, B, C and T* are the usual Landau-De Gennes 
parameters. One may note the explicit dependence of 
Eq. H^9(l on the short-length cutoff a and be concerned 
about the sustainability of the continuum limit a — > 0. 
In fact, because g oc a 3 and (3k ~ 1/a, the powers of a 
cancel and the free energy does not depend on this cutoff 
parameter, which is a reassuring test of consistency of 
our theory. 

The easiest way to illustrate the consequence of the 
free energy renormalization is to plot the equilibrium 
values of the order parameter Q. The values chosen 
for the Landau phenomenological constants are: Ao ~ 
5.0 x 10 3 Jm^Jf' 1 , B ~ 3.3 x 10 5 JnT 3 and C ~ 
1.0 x I0 6 Jm~ 3 . A detailed description on how they are 
obtained by analyzing experimental data of Ref. [22j is 
given in the Appendix. We consider these values to be 
indicative only since different nematic materials will cer- 
tainly have large variations in these values. Plotting 
the equilibrium order parameter against reduced tem- 
perature, Fig. |3 we can compare Q(T) for the disorder 
strength increasing from g = in (a) up to g — 10~ 4 
in (d). As mentioned in Section IiTtI a typical polydo- 
main elastomer should have g ~ 4 x 10~ 3 and therefore 
the model predicts a supercritical behavior in agreement 
with many experiments. Note that in contrast to this re- 
sult, the quenched orientational disorder was shown not 
to alter the continuous nematic phase transition in thin 
films, whose director is confined in the AF-plane [38| . 

The inclusion of disorder has a profound effect on the 
phase behavior of 3D systems whose pure versions un- 
dergo a first order phase transition. The discontinuous 
jump of the order parameter at the nematic transition 
becomes smaller as the strength of the disorder g in- 
creases and eventually disappears altogether above a crit- 
ical value making the phase transition continuous. The 



Q 




I ■ L — , — , — , — , — , — , — , — , — > 

1 1.02 1.04 1.06 



FIG. 2: The equilibrium order parameter Q as a function of 
reduced temperature t = T/T* for a range of different disor- 
der strength g. (a) is the first order transition for a system 
with no disorder (b) is a subcritical system (c) is a critical 
system (d) is a supercritical system. As the disorder strength 
increases the discontinuous jump decreases and eventually dis- 
appears. 




FIG. 3: Free energy against order parameter plots for a sub- 
critical system for a range of temperatures, Ti < T2 < T3 < 
T4. The doted and dashed lines show respectively the Lan- 
dau and disorder energies for T = T3. They illustrate that 
the high-Q minimum is a product of the Landau energy exclu- 
sively, whereas the low-Q minimum arises from competition 
of the A Q 2 Landau term and the disorder energy. 



change in behavior is explained by the simple fact that 
the energy terms arising from disorder scale as negative 
powers of Q. As a direct result, the energy of the system 
increases as the transition approaches and the zero order 
parameter phase is never reached. To see exactly how 
this happens consider plots of F(Q) of a subcritical and 
a supercritical system. 

For a subcritical system, where disorder is weak, the 
jump in Q(T) is still present, albeit smaller than in the 
original transition of the pure system. Because param- 
eter g is small, the disorder part of the energy becomes 
significant only for small Q, where it diverges. The ap- 
pearance of the jump has the same origins as in the clas- 
sical Landau-De Gennes theory. Fig. [3] shows plots of 
the energy density against order parameter for four dif- 
ferent temperatures around T/yj. At the lowest tempera- 
ture (Ti) the single minimum determines the equilibrium 
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value of Q. As the temperature increases the high-Q 
minimum moves to a slightly smaller value of Q and - 
more importantly- another "low-Q" minimum appears as 
a result of disorder. At T2 the high-Q is still the global 
minimum, but at the critical temperature (T3) the two 
minima have the same energy value. This means that 
two distinct phases, one with Q ~ 0.23 and the other 
Q ~ 0.05, coexist. Once this temperature is passed (T4) 
the low-Q minimum determines the system's order pa- 
rameter. The crucial difference between these plots and 
the classical Landau-De Gennes theory is that the low-Q 
minimum in the latter is always placed at Q — 0. Since 
the disorder energy diverges at zero Q, this minimum is 
pushed at positive values of Q in systems with quenched 
disorder. 

The dotted and dashed lines show the Landau and 
the disorder energy for the same temperature T3, respec- 
tively. Around the high-Q minimum the dotted line has 
the same shape as the actual energy; apart from a con- 
stant shift to lower energy they are exactly equal. There- 
fore this minimum is a result of the competition of the 
—B Q 3 and C Q energy terms of the Landau expansion, 
which dominate at large Q. The position of the low- 
Q minimum is influenced by disorder. This minimum 
is a balance of the divergent Fd term and A Q 2 . For 
temperatures well above Tjvi the minimum is located at 
very small order parameter and the only relevant terms 
of Eq. dSBJ) are DQ~ 4 and AQ 2 . 

In a supercritical system disorder is stronger (large g) 
and therefore its effect on the energy is more prominent. 
As a result the effect of Fd is relevant for all the values 
of Q, not just in the small order parameter region as in 
the previous case. Fig.^Jshows the relevant energy plots. 
Crucially there is only one minimum at any given tem- 
perature. Its position shifts to smaller order parameter 
as temperature increases, but the phase transition is con- 
tinuous. In comparison with Fig. [3J we can say that the 
low-Q minimum has "broadened" and "absorbed" the 
high-Q minimum. 

Another difference with the small g case is that, be- 
cause Fd has larger magnitude, the low order A Q 2 term 
of the Landau expansion does not influence the position 
of the minimum. The dotted lines in Fig. 0] are drawn for 
the same temperature as in the previous case of Fig. |3 
but now this temperature is labeled as T\ . In this super- 
critical system the energy minimum occurs at approxi- 
mately the same value of order parameter as the high- 
Q minimum of the Landau DeGennes expansion (thin- 
doted line). In a pure system this would be a metastable 
state because there would exist a global minimum at 
Q = 0. 

The natural question to ask next is what is the critical 
point at which the jump in Q(T) disappears completely. 
At this point the free energy must have the merging of 
all its minima and maxima. Solving for the first, second 
and third derivatives being zero provides three equations 
for the unknown critical parameters g c , T c and Q c and 
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FIG. 4: Energy against order parameter plots of a supercrit- 
ical system for a range of temperatures, Ti < T2 < T3, 
all larger than Tnj. The dotted and dashed lines show the 
Landau and disorder energies for Ti, respectively. Since the 
disorder energy diverges the AQ 2 Landau term is no longer 
significant. 



the critical point is then given by 

(9c Qc) 



(2^ 2 7 11 B 12 { K [3af_ 7B_\ 



(40) 



It is worth comparing this with the classical case of an 
applied field, which adds a —fQ term in the Landau 
expansion, where the critical Q is slightly smaller and 
equal to B/3C. 

In the analysis of the replica symmetry stability in 
Section |W] we found that the above description of ne- 
matic systems breaks down for small order parameter 
when Q < Q s t a b, see Eq. (|57|) . The threshold value of 
Qstab is proportional to the cube of the cut-off length 
a. In order to be consistent with the line adopted in 
previous calculations, the continuum limit of space was 
employed and a was taken to zero. Hence Q s tab was also 
taken to be zero. However, in reality the length scale a 
is in fact the nematic coherence length, because below 
this size we cannot write Frank elasticity and there is no 
meaning to order parameter Q, or director n. Keeping 
a non-zero, and assigning it the value of 5nm which is 
usually associated with liquid crystals, makes Q s tab also 
finite. Nevertheless substituting the value of g ~ 3 x 10 -3 
and k ~ ksT/a gives Q s tab — 2 x 10 -4 . This is an ex- 
tremely small value for the order parameter and, as we 
see from Fig|3 it would only be acquired at temperatures 
well above the experimental range. Hence, even if a =/= 0, 
the window of replica symmetry stability is wide enough 
to describe realistic nematic systems. 



VI. SUMMARY 

This paper examines how the inclusion of randomly 
quenched orientational disorder leads to the rounding of 
the nematic-isotropic phase transition in three dimen- 
sions. The coupling between impurities and the local or- 
der parameter pins some mesogenic molecules and does 
not allow the sample to have a uniform director field n(r). 
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After quenched disorder has been averaged over using 
the replica method, a replica-symmetric auxiliary field is 
used to obtain the free energy arising from disorder. The 
disorder energy adds a oc Q~ 4 term to the Landau-De 
Gennes expansion which diverges for diminishing order 
parameter Q. As a result the isotropic phase is never 
reached and, for sufficiently strong disorder, the phase 
transition becomes continuous. This is in accordance 
with many experiments on nematic elastomers that also 
show a smooth transition rather than a discontinuous as 
predicted by the classical Landau-De Gennes theory. An 
earlier study on XY nematics, whose director is confined 
to a plane, showed that the quenched disorder does not 
affect their continuous phase behavior |S8j . 

A stability analysis shows that the replica-symmetric 
solution we have employed does fail at small values of 
the order parameter, at Q < Q s tab- However, for realis- 
tic values of physical parameters we estimate the order 
of Qstab ~ 1CP 4 in nematic elastomers. A supercritical 
system acquires such low values of Q at temperatures 
well above the elastomers melting point and therefore 
the window of replica symmetry stability adequately de- 
scribes the 3D nematic system. 
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APPENDIX 

In order to determine the three phenomcnological pa- 
rameters Aq, B and C of the Landau-De Gennes the- 
ory, defined in Eq. (|38|) , we need three independent mea- 
surements. These are provided by a NMR experiment 
measuring the order parameter as a function of tempera- 
ture of a polymer melt and its corresponding crosslinkcd 
network (a polydomain nematic elastomer) |22|. As ex- 
pected from the Landau-De Gennes theory, the melt 
shows a first order transition with the discontinuous jump 
in order parameter being approximately AQ ~ 0.22, 
which is smaller than 0.4, the value usually associated 
with ordinary liquid crystals |3^|. The theoretical pre- 
diction gives this jump to be equal to AQ = 2B/3C 
|14| and substituting the experimental measurement we 
find B ~ 0.33C. The second measurement is the width 
of the temperature hysteresis which is AT ~ 5K in the 



polymer melt. Theory predicts AT = 2B 2 /9A C. Com- 
bining this with B ~ C/3 we find B = 67 Aq IK and 
C = 200 Ao IK. To get a value of Aq a third measure- 
ment is required. 

A striking difference between the two Q{T) plots, for a 
nematic polymer melt and its crosslinked elastomer ver- 
sion, is that crosslinking has reduced the overall order 
at temperatures below T/yj. For elastomers crosslinked 
in the isotropic phase, the energy addition arising from 
nematic rubber elasticity adds a fourth order term in the 
Landau expansion (Ipf : 



-/*a Q 



where \i is the rubber modulus and a accounts for the 
microscopic details of an elastomer. For a freely joined 
polymer a = 3 [TlT ] , but a side-chained polymers have a 
ranging between —0.5 and |4jJ. Let us take an inter- 
mediate case where a = 1. When this term is added to 
the Landau-De Gennes expansion, the fourth-order coef- 
ficient (which is C/4 for the polymer melt) now becomes 
larger x(l + 3/ia 4 /C). Hence the transition tempera- 
ture, given by T NI = T* + 2B 2 /3A C 2 , decreases since 
the renormalized C increases. The shift in this transi- 
tion temperature between the melt and the correspond- 
ing elastomer is: 



AT NI = - 



2BV 
3A C 2 



(A.l) 



and it provides the third relation that allows to determine 
Aq. An estimate of ATni is possible in Ref j^: Tni is 
easily identified in the polymer melt, but it is not clear 
what it means in the disordered nematic elastomer with 
a continuous transition. It can be loosely defined as the 
temperature where Q = AQ = 0.22. This then makes 
AT/vi ~ 15° C. Substituting this back to Eq. ifO) with 
a = 1, we obtain Aq ~ 200 K ' ^ typical nematic elas- 
tomer has elastic modulus of the order of 10 6 Pa. Putting 
everything together, the phenomenological constants of a 
nematic elastomer are: 



Aq 



5.0 x 10 3 Jm" 3 ^ -1 , 



B 



3.3 x 10 5 Jm" 3 



and C 



1.0 x 10 6 Jm -3 . 



These values are crude estimations, only given here to 
illustrate the effect of disorder in our model. However, 
it is comforting that, although obtained from a different 
set of experimental measurements, these values are quite 
close to the ones reported in literature pd| . 
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